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1. Introduction 

With the advent of particle accelerators in the few GeV energy region, 
heavy meson production in hadronic collisions has attracted increasing at- 
tention in the past few years. In particular, heavy meson production in 
nucleon-nucleon (NN) collisions is of special interest because it allows to 
investigate in a simple system the short range hadron dynamics for which, so 
far, we have a very limited knowledge. Due to the large momentum trans- 
fer between the initial and final states, these reactions necessarily probe 
the dynamics at short distances. In this context, apart from the intrinsic 
interest associated with the particular meson produced, the production of 
pseudoscalar and vector mesons can be used as a tool to probe the short dis- 
tance dynamics systematically. Table p] illustrates this point: it shows the 
momentum transfer and the corresponding distance probed by producing 
mesons of different masses at the respective threshold energy. At threshold, 
the momentum transfer is given by q = (m^mM+'ni 2 M /A) 1 /' 2 , where mjy and 
ttim denote the nucleon and meson mass, respectively. As one can see, the 
distance probed in these reactions ranges from 0.53/m for pion production 
to 0.18/m for meson production. 



Table 1. Momentum transfer q and the corresponding distance r probed by the 
NN — > NNM reaction at the threshold energy for different particles M produced. 



particle 


mass (MeV) 


q (fm x ) 


r (fm) 


7 





0.0 


oo 


7T 


140 


1.9 


0.53 


V 


550 


3.9 


0.26 


p,u 


780 


4.8 


0.21 


rf 


960 


5.4 


0.19 


4> 


1020 


5.6 


0.18 



The theory of heavy meson production is still in its early stage of de- 
velopment. As we have seen above, heavy meson production reactions can 
probe quite short distances - down to less than 0.2/m. These short dis- 
tances correspond to the region of confinement where the relevant degrees 
of freedom are the constituent quarks and gluon flux tubes. Thus, the ap- 
propriate approach to describe these short range dynamics might be the 
constituent quark models rather than the hadronic models. However, such 
an approach still remain to be developed. On the other hand, the transition 
region from the hadronic to constituent quark degrees of freedom does not 
have a well-defined boundary and, consequently, it is of special interest to 
see how far down in distance one can "push" the hadronic models. 
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In principle, effective field theories can be used to describe near-threshold 
particle production processes in terms of hadronic degrees of freedom. How- 
ever, although such an approach (Chiral Perturbation Theory) might be 
appropriate for describing the production of pions ||] in spite of the rela- 
tively large expansion parameter, Q = y / m^/m/v ~ 1/3, for the description 
of heavy meson production there is, a priori, no obviously preferred ap- 
proach. The majority of existing calculations of heavy meson production in 
NN collisions are based on meson exchange models of strong interactions 
[^, [|, |3|, ^ q, 0, g, 1C , 11, 12] with a few exceptions [13|. The price one pays 
for insisting on such models is that their predictions become more and more 
sensitive to the short range part of the model, usually parametrized in terms 
of the form factors at the hadronic vertices involved. The success of such 
models should be measured in terms of their capability to correlate as many 
independent processes as possible in a consistent manner. We mention that 
although the approach used in Ref.[12] is based on meson exchange models, 
it differs from conventional meson exchange models in a number of aspects, 
such as the absence of the form factors at the hadronic vertices, etc. 



2. Meson Exchange Models 

As mentioned in the previous section, the majority of the existing cal- 
culations of meson production in iViV collisions are based on conventional 
meson exchange models. Even within such models, the description of these 
reactions is not a simple task in principle, for the final state is a three-body 
state and, consequently, one needs to solve the three-body Faddeev equa- 
tion. Of course, a complete three-body calculation of this reaction is at 
present not available [pH| . How can one then start to describe these reac- 
tions? In order to gain some insight to this question, let us examine some 
of the features exhibited by the production cross sections near threshold. In 
reactions like AW bremsstrahlung, where a (massless) photon is produced, 
the measured cross section varies with the inverse of the photon energy w 7 
near threshold. This feature of the cross section is expressed by the so-called 



soft-photon theorem [15] which gives 



a = — + B, (1) 

where A and B are constants containing only the on-shell information of 
the NN interaction (or, in other words, the asymptotic behavior of the AA 
wave function). This result is not surprising at all if we recall from Tabic 
|l| that, close to threshold, the photon production reaction in AA collisions 
probes only the asymptotic behavior of the NN wave function. Eq.(|l]) is, 
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therefore, regarded as a model-independent result and, as such, it should be 
reproduced by any model describing the NN bremsstrahlung reaction. 

For production of heavy mesons, we do not have a model independent 
result like the low energy theorem for NN bremsstrahlung. However, as 
early as 1952, Watson |l6[ showed that in heavy particle production reac- 
tions in which strongly (attractively) interacting particles are present, the 
energy dependence of the cross sections should be dictated by the energy 
dependence of the interaction between those strongly interacting particles 
and the available phase space. In particular, for meson production in NN 
collisions, the energy dependence of the total cross section should be given 
by the energy dependence of the on-shell NN final state interaction (FSI), 
T(p',p'), plus the phase space 

a(E) oc / dp(E,p')\T(p',p')\ 2 




where the integration is over the available phase space, p(E,p'), with p' 
denoting the relative momentum of the two interacting nucleons in the final 
state. 5(p') denotes the corresponding NN phase shift. Watson's result is 
based on the observation that the massive particle production is a short 
range process and, as such, the primary production amplitude of such a 
particle should have a weak energy dependence. Note that the production 
cross section cr(E) cannot be expressed in a model-independent way, for its 
absolute value depends on the short range part of the interaction. We shall 
elaborate more on Watson's result later. For the moment, we mention that 
all the recently measured meson production cross sections in NN collisions 
near threshold follow the energy dependence given by Eq.(Q) with the ex- 
ception of T) meson production, where one sees a relatively small deviation 
at energies close to threshold. This deviation is commonly attributed to the 
r)N FSI. 

The above consideration indicates that any model of heavy meson pro- 
duction reaction in NN collisions should have built in the feature given 
by Watson's result, Eq.(^). As we shall show in the following, this can be 
achieved within a Distorted Wave Born Approximation. Here we follow a 
diagrammatic approach to derive the total amplitude. We start by consid- 
ering the meson-nucleon {MN) and NN interactions as the building blocks 
for constructing the total amplitude describing the NN — > NNM reaction. 
We then consider all possible combinations of these building blocks in a 
topologically distinct way, with two nucleons in the initial state and two 
nucleons plus a meson in the final state. In this process of constructing 
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the total amplitude, care must be taken in order to avoid diagrams that 
lead to double counting. Specifically, these are the diagrams that lead to 
the mass and vertex renormalizations, since we choose to use the physical 
masses and coupling constants. The resulting amplitude constructed in this 
way is displayed in Fig. |l|. The ellipsis indicates those diagrams that are 




(1< — >2) + o o o 

Fig. 1. Meson production amplitude obtained in a diagrammatic approach to the 
NN — > NNM reaction considering the MN and NN T-matrices as the basic 
building blocks. Tmn stands for the MN T-matrix. FSI and ISI stand for the final 
and initial state NN T-matrices, respectively. The first diagram on the r.h.s. is 
referred to as the production current J which enters in other diagrams as can be 
noted. 

more involved numerically (including, in particular, the MN FSI, which 
otherwise would be generated by solving the three-body Faddeev equation). 
So far there are very few attempts to account for them || |I7|] . Therefore, 
neglecting those diagrams, the total amplitude is given by the diagrams 
displayed explicitly in Fig. |l] and reads 

M = (1 + T\- )] G { f y ) J(l + a[ +) T/ +) ) , (3) 

where Tuj\ denotes the AW T-matrix interaction in the initial(i)/final(/) 
state and, Gnj\, the corresponding two nucleon propagator. The super- 
script ± in Tuj\ as well as in Gr(i,/) indicates the boundary condition ((— ) 
for incoming and (+) for outgoing waves). The production current is de- 
noted by J, which is nothing other than the MN T-matrix, Tmn, with 
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one of the meson legs attached to a nucleon (first diagram on the r.h.s. in 
Fig. H). Eq.(H|) is the basic formula on which the large majority of existing 
calculations are based. 

We now wish to exhibit the essential features of the meson production 
reaction contained in Eq.(^). In particular, we would like to make close 
contact between the amplitude M and Watson's result given by Eq.(^). To 
this end, we use the relation 

where a = i, /, to express M as 

M = - iKf T f (p\p')[l + i-^V f (E,p')}} J(E,p') 

x \l-i Ki Ti(p, P )[l + i— Vi{E,p')] \ , (5) 
I ap J 

where Ti(p,p) denotes the on-shell NN T-matrix with the relative momen- 
tum p of the interacting two nucleons in the initial state i and Tf(p',p') 
denotes the on-shell NN T-matrix with the relative momentum p' of the 
interacting two nucleons in the final state /. Here the superscript (+) of 
Tuf\ has been omitted. Ki = irpe(p)/2 and kj = irp'e(p')/2 are the phase 
space densities of the two nucleons in the initial and final states, respec- 



tively. e(q) = y q 2 + m N and a is a constant (scattering length) introduced 

for convenience. Note that in the above equation, only those arguments 
of the quantities relevant to the discussion related to Watson's result are 
explicitly displayed. E = E{. The function V% is given by 

V iM = f2V r^J^L, (6a, 







f(F h] _ Tj(p,k)J(E,k) _ Kj{p,k)J{E,k) 



Similarly, 



T f (p',k)A(E,k) K f (p>,k)A(E,k) 
Jf{ ' ' T f (p',p>)A(E,p>) K f {p>,p>)A{E,p') ' 1 j 

with -f^(ij) denoting the NN K-matrix and A = (1 + GjTj) J. The functions 
V^j^Ejp') summarize all the off-shell effects of the NN interaction and 
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production current. As such, they are unmeasurable and model-dependent 
quantities. For production of heavy mesons near threshold, the off-shell 
NN interaction required in Vf(E,p') is at very low energy. In calculations 
based on conventional meson exchange models this function is very large and 
cannot be neglected. In contrast, the function (1 / ap)Vi{E , p') is relatively 
small. This is because for heavy meson production the energy of the two 
nucleons in the initial state has to be large enough to produce the meson in 
the final state. For example, for r\ meson production, the incident energy of 
the beam nucleon corresponding to the threshold energy is about 1.25 GeV. 
At such high energies, the on-shell NN interaction has a rather weak energy 
dependence. Note that the inverse scattering theory tells us that, for a 
local and energy-independent NN potential, the off-shell behavior of the 
NN amplitude is completely determined if one knows the corresponding 
on-shell amplitude in the entire energy domain. At least in the case of 
meson exchange models, the weak energy dependence of the on-shell NN 
interaction implies also a flat off-shell behavior and leads to a small value 
of (l/ap)Vi(E,p'). Therefore, at least for the discussion of the essential 
features of the cross section near threshold, the function (l/ap)Vi(E,p') 
may be neglected. 

Using then the relation between the on-shell T-matrix and the phase 
shift 5(q) and inelasticity rj(q) 

«(q)T(q,q) = \ (ry^e 2 ^ - l) , (8) 
Eq.(|5|) can be reduced to 

M <* |- e W (^^) [l + V f (E,p')]}j(E,p') 

x {^(P)e^ (p) + l)} , (9) 

where the inelasticity in the final AW state is set to unity, Tjf(p') = 1. 
The above equation is the desired result. It exhibits the essential features 
of the meson production reaction in NN collisions. First, it shows that 
the relevant energy dependence of the total amplitude near threshold is 
indeed determined by the strong energy dependence of the on-shell NN 
FSI, proportional to sm{5f{p'))/ap' . This, when combined with the phase 
space factor, determines the energy dependence of the cross section as given 
by Eq.(|2|). Note that for heavy meson production the production current 
J(E,p'), as well as Vf(E,p'), should be weakly energy-dependent, for they 
summarize all the short range dynamics of the reaction. Thus, none of the 
terms in Eq.(^), apart from the on-shell NN FSI, should introduce any sig- 
nificant energy dependence near threshold; they just amount to a constant. 
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Therefore, the total cross section data near threshold determine just this 
constant. Second, the major effect of the NN initial state interaction (ISI) 
(the term in the last curly brackets) is to reduce the cross section by a factor 
of 



A -[, /((/ ,„.'-' MW + l 

= m ( P ) cos 2 (^( P )) + i[i - Vi ( P )} 2 < i[i + ^(p)] 2 . (10) 

As mentioned before, the majority of the existing calculations of meson 
production in NN collisions are based on Eq.(||). The differences among 
them reside in how the production current J is modeled, as well as in the 
different treatments of both the NN FSI and ISI. Let's first concentrate 
on the NN interactions. It is clear that the description of meson produc- 
tion processes in NN collisions based on Eq.(H|) requires the half-off-shell 
NN ISI and FSI. This has been exhibited explicitly in Eq.([|) through the 
functions Vuj\ given by Eqs.(||,|7]). Although the on-shell NN interaction 
can be determined from the NN elastic scattering experiments, the off-shell 
behavior of it can only be provided by a given model of the NN interaction. 
For energies below pion threshold, there exist a number of accurate meson 
exchange models [18, [0^, 2C] which can provide the half-off-shell extension 



of the NN interaction. As previously noted, these so-called realistic NN 
interactions, based on meson exchange models, yield a rather large value 
for the function Vf(E,p') that, consequently, cannot be neglected. In this 
connection, for production of heavy mesons, the predicted total cross sec- 
tions can easily differ by a factor of two or more due to the different off-shell 
behavior of these realistic NN FSIs. This indicates that a consistent treat- 
ment of the NN FSI and the production current J is required. Maintaining 
the consistency between the NN FSI and the production current, however, 
is not a trivial task. While in models where the underlying meson exchange 
structure is known this consistency can, in principle, be maintained, in the 
case of purely phenomenological models, such as the parametrized version of 
the Paris NN interaction |ll|, such a consistency is impossible to achieve - 
even in principle. In the excess energy region below Q = 100 MeV however, 
the introduced difference in the predicted total cross section due to off-shell 
differences of these realistic NN interactions is practically a constant. It 
should also be mentioned that the procedure of just evaluating J in the on- 
shell tree level approximation and simply multiplying it by the on-shell NN 
FSI, as has been done by many authors (see references quoted in plf ), is not 
acceptable for obtaining quantitative predictions. As it can be seen from 
Eq.(|9|), the strength of the amplitude M depends on the function Vf(E,p'). 
Using just the on-shell T-matrix as the FSI instead of the full half-off-shell 
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T-matrix means setting the function T , f(E,p') to zero. Such a procedure 
simply lacks a consistent regularization scheme between the NN interaction 



and the production current J [2jJ. 



One of the major limitations in the current models of heavy meson 
productions in NN collisions is the lack of a reliable model for the NN 
ISI. As mentioned before, the energy of the two interacting nucleons in the 
initial state must be large enough to produce a meson. For example, for 
the r\ meson this means nucleon incident energies of at least 1.25 GeV . For 
heavier mesons, like rf and (f), the corresponding threshold incident energy is 
well above 2 GeV, where no reliable model exists to provide the half-off-shell 
NN T-matrix that is required in the evaluation of the function Vi(E,p') in 
Eq.@. Note that, although this quantity is expected to be small and may be 
neglected for estimates of cross sections, for more quantitative predictions it 
should be taken into account. In particular, predictions of spin observables 
are expected to be sensitive to this function. The situation with the NN ISI 
is even worse in the case of heavy meson production in pn collisions. There, 
we lack even the on-shell NN interaction. While for total isospin T = 1 
states rather reliable phase shift analyses exist up to 3 GeV incident energy, 



for T = states the reliability is limited to 1.3 GeV [22, 23|. This situation 
imposes a severe limit in all existing models of production of mesons heavier 
than the rj meson in pn collisions. 

There are basically three different approaches in literature in modeling 
the production current J based on meson exchange models. One is a mi- 
croscopic model of MN — > M'N' reactions whose data are reproduced by 
the model ||, pj] ]. The off-shell behavior of the MN T-matrices that are 
required in the construction of J is then provided by the model. Another 
approach is a "pseudo" empirical model in which the production current is 
constructed using the on-shell amplitudes of the MN — ► M'N' as well as 
jN — > M'N' reactions extracted directly from the available data [l(| . In 
the latter reaction the VMD is used to convert its on-shell amplitude to an 
MN — ► M'N' amplitude (M = p,u). The off-shell behavior of the corre- 
sponding amplitudes necessary to construct J is then an assumption in this 
approach. The third approach is to split the MN T-matrix into the pole 
(T p ) and non-pole (T NP ) parts and consider the non-pole part in the Born 
approximation only in order to construct the production current J || 
As has been shown elsewhere ^5|, the MN T-matrix can be split into the 
pole and non-pole parts according to 

T = T p + T NP , (11) 

where 

T P = f M NB9BfMNB , (12) 
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with /mnb and gs denoting the physical meson-nucleon-baryon (MNB) 
vertex and baryon propagator, respectively. The summation runs over the 
relevant baryons B. The non-pole part of the T-matrix is given by 

T NP = yNP + yNP GT NP (13) 

where V = V — V p , with V p denoting the pole part of the MN potential 
V. V p is given by the equation analogous to Eq.(p"^) with the renormalized 
vertices and propagators replaced by the corresponding bare vertices and 
propagators. This third approach then leads to a production current which 
is obtained by approximating the full MN T-matrix as T = T p + V NP . 



3. Our Model 

In our model of NN — > NNM, the reaction amplitude M is calculated 
using Eq.@. It is based on a relativistic meson exchange model of hadronic 
interactions. The results we shall show in the next section are obtained by 
using the NN interaction developed by the Bonn group [ |18| . This interac- 
tion is obtained by solving a three dimensionally reduced (Blankenbecler- 
Sugar) version of the Bethe-Salpeter equation. The loop integrals in Eq.(^|) 
are calculated consistently with the three dimensional reduction used in con- 
structing the NN interaction, i.e., the two-nucleon propagators Guj\ are 
taken to be the Blankenbecler-Sugar propagator. The production current J 
is modeled according to the third approach discussed in the previous section. 
It consists of the nucleonic, resonance and mesonic currents as displayed di- 
agrammatically in Fig. |2[ Note that they are all Feynman diagrams and, 
as such, they include both the positive- and negative-energy propagation of 
the intermediate particles. 

All the parameters of our model for the production current (coupling 
constants and form factors) are confined to the hadronic vertices Tmnb, 
^mnb and Tmmm indicated in Fig. ^ (B = N,N*). Below, we discuss 
briefly these parameters for each vertex. For more details, see Refs.|7], B. 

Fmnb'- At the meson production vertex Tmnn, the coupling constant is taken 
consistently with that used in the construction of the NN interaction. 
In addition, this vertex contains an extra form factor (Fj^(p 2 )) to 
account for the off-shellness of the intermediate nucleon, which is far 
off-shell if a heavy meson is produced. It is given by 
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(1 <— > 2) 



Fig. 2. Model for the meson production current J. The upper two diagrams 
are called nucleonic or resonance current depending on the intermediate baryon 
B being a nucleon (N) or a nucleon resonance (N*). The lower diagram is the 
mesonic current. M, M' = tt, rj, p, lu, a, a Q . 



with Atv = 1.2 GeV. p and tun stand for the four-momentum and 
mass of the intermediate nucleon, respectively. 

The coupling constant in the vertex Tmnn* is extracted from the 
measured decay width of the resonance into a meson and a nucleon, 



N* — > M+N, whenever available [26]. For the vertex involving a vec- 



tor meson, the coupling constant can be extracted from the measured 
radiative decay using the VMD. The sign of the coupling constant 
is chosen in accordance with the relevant photo-production reaction 



analysis p7l 28 1. Tmnn* is a ls° multiplied by the form factor given 



by Eq.(ffl), with tun replaced by m^v*. 

^mnb- The vertex Tmnn is taken consistently with that used in the con- 
struction of the NN interaction. For Tmnn involving the off-shell 
nucleon that produces the meson, it is also multiplied by the form 
factor Fn(p 2 ) given by Eq.(|l4|). The vertex Tmnn* is the same as 
Tmnn* j except that it contains an extra form factor due to the off- 
shell meson which is taken to be the same as that of the corresponding 
vertex Tmnn- 



Tmmm- The coupling constants in the three-meson vertices (Tmmm) are ex- 
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tracted from both the strong and radiative (measured) decay widths 



[26] in combination with SU(3) symmetry, plus the OZI rule |pS|| . The 
latter relates the basic SU(3) octet and singlet coupling constants. 
Fmmm includes the form factor 




with Am = Am' = 1.45 GeV. q and q' denote the four momenta of 
the two mesons with masses tum and mj,/', respectively, that fuse to 
produce the third meson. The parameter x(= 0, 1) ensures the proper 
normalization according to the normalization point q' 2 = 0, mL at 
which the corresponding coupling constant is extracted . 

There is a number of features in the model described above which are per- 
haps worth mentioning here. This model is suited for a systematic analysis 
of the production of different mesons within the same model due to the 
simplicity of modeling the production current. Also, the model is especially 
suited for studying the role of different reaction mechanisms that produce 
a meson. The consistency between the NN interaction and the production 
current can be easily maintained when one knows the underlying meson 
exchange structure of the NN interaction used (note that the way the pa- 
rameters of our model are fixed ensures the consistency between the NN 
interaction and the production current). As mentioned before, this is critical 
if one wishes to achieve quantitative predictions, especially for production 
of heavy mesons. 



4. Some selected results 



In this section we shall present some selected results on the vector and 
pseudoscalar meson productions based on the model described in the pre- 
vious section. As an example of vector meson production, we consider <p 
meson production in pp collisions. For the pseudoscalar meson production, 
we discuss the results for rj production in both the pp and pn collisions. 



4-1- pp — > pp<j) 



The study of this reaction is of particular importance in connection 
with the questions related to the amount of hidden strangeness in the 
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nucleonjEO, [31] ] . In the context of (ft meson production processes one ex- 



pects [31] that a large amount of hidden strangeness in the nucleon would 



manifest itself in reaction cross sections that significantly exceed the values 



estimated from the OZI rule [29]. This phenomenological rule states that 
reactions involving disconnected quark lines are forbidden. In the naive 
quark model the nucleon has no ss content, whereas the (ft meson is an ide- 
ally mixed pure ss state. Thus, in this case, the OZI rule implies vanishing 
<j)NN coupling and, accordingly, a negligibly small production of (ft mesons 
from nucleons by non-strange hadronic probes. In practice there is a slight 
deviation from ideal mixing of the vector mesons, which means that the (ft 
meson has a small uu + dd component. Thus, even if the OZI rule is strictly 
enforced, there is a non-zero coupling of the (ft to the nucleon, although 
the coupling is very small. Its value can be used to calculate lower limits 
for corresponding cross sections. For example, under kinematic conditions 
chosen to cancel out phase space effects, one expects cross section ratios of 
reactions involving the production of a (ft- and an u> meson, respectively, to 
be 

a(A + B —> <j>X) 2 

R = (A,n V\ ~ tan ay ) ' ( 16 

a [A + B — > ujX) 

where A, B and X are systems that do not contain strange quarks, ay = 
By — By (ideal) is t ne deviation from the ideal u — (ft mixing angle. This 
result arises from simply equating the cross section ratio to the square of 
the ratio of the relevant coupling constants at the (ft and to production 
vertices. According to the OZI rule plus SU(3) symmetry, these couplings 
are proportional to sin(ay) and cos(ay), respectively. With the value ay = 



3.7° [26] one gets a rather small ratio of R = 4.2 x 10 . The data presented 



by the DISTO collaboration 32, 33 1 in pp collisions indicate that this ratio, 
after correcting for phase space effects, is about eight times larger than the 
above OZI estimate. 

In principle, the (ft meson production cross section can be used for a 
direct determination of the (ftNN coupling strength. Any appreciable (ftNN 
coupling in excess of the value given by the OZI rule (g^NN — — (0.60±0.15)) 
[0] could be interpreted as evidence for hidden strangeness in the nucleon. 
Of course, there is also an alternative picture: one in which the coupling of 
the (ft meson to the nucleon does not occur via possible ss components in 
the nucleon, but via intermediate states with strangeness. Specifically, this 
means that the (ft meson couples to the nucleon via virtual AK, T,K, etc. 
states. Corresponding model calculations |34], ^] have shown, however, that 
such processes give rise to (effective) (ftNN coupling constants comparable 
to the OZI values and therefore should not play a role in drawing conclusions 
concerning hidden strangeness in the nucleon. 



14 



CRACOW PRINTED ON FEBRUARY 8, 2008 



Details of our calculation may be found in Ref. J5J . The work of Ref . pi] is 
a combined analysis of the u and (j) meson production in pp collisions in order 
to reduce the number of free parameters in the model. Here we only report 
on the essential results. For the production current J we consider only 
the nucleonic and mesonic currents. The latter consists of vpir exchange 
current (v = u),(f>). Other mesonic currents have been investigated in a 
systematic way and found to be very small. The nucleon resonance current 
is not considered because, presently, there is no established resonance that 
decays into a vector meson and a nucleon (see Ref. |36| in this connection) . 
It is found that the (j) meson is produced significantly from the (ppir exchange 
current in pp collisions. This means that one must be able to disentangle this 
current from the nucleonic current (where the (j) meson is emitted directly 
from the nucleon) if one wishes to extract the (j)NN coupling strength from 
this reaction. As pointed out in Ref.H, this can be done by looking at the 
angular distribution of the emitted meson. Because the mesonic current 
contribution yields an isotropic angular distribution, whereas the nucleonic 
current contribution exhibits a cos 2 (#) dependence, the two currents can be 
disentangled uniquely from the angular distribution data. Fig. || illustrates 
this point for the pp — > ppuj reaction at a proton beam energy of 2.2 GeV . 
The upper panel illustrates the situation when the nucleonic current (dashed 
curve) is smaller than the mesonic current (dash-dotted curve) resulting in 
a nearly flat angular distribution (solid curve). Note that the interference 
between the two currents is destructive, which is a direct consequence of 
the signs of the relevant coupling constants determined in our model as 
discussed in the previous section. The lower panel shows the situation when 
the nucleonic current is larger than the mesonic current. Here the angular 
distribution exhibits a pronounced angular dependence. In both scenarios 
the total cross sections have been kept to be about the same. 

The earlier data of the 4> meson angular distribution by the DISTO col- 
laboration (3^] at a proton beam energy of 2.85 GeV is shown in Fig. ||] 
together with our result. The absolute normalization of this data has been 
determined as described in Ref. |7j . Although the data have large uncertain- 
ties, the observed angular distribution is rather flat. Note that the angular 
distribution should be symmetric about 9 = 90°, due to the identity of the 
two protons. The solid curve is one of our calculations fitted to the data. 
The dash-dotted and dashed curves are the corresponding mesonic and nu- 
cleonic current contributions, respectively. As one can see, the data require 
a very small contribution from the nucleonic current. The value of the (pNN 
coupling constant thus extracted is in the range g<f>NN ~ —(0.2 — 0.9). More 
data for both the <p an d oj meson at low excess energies are necessary in 
order to determine better the parameters of the model, and thus reduce the 
uncertainties in the extracted value of g^NN- I n this connection, current ex- 
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Fig. 3. Angular distribution of the emitted uj meson in the cm. frame of the total 
system for two possible scenarios at a proton incident energy of 2.2 GeV . The 
upper panel corresponds to the case when the nucleonic current is smaller than the 
mesonic current. The lower panel corresponds to the case when the nucleonic cur- 
rent is larger than the mesonic current. The dashed curves represent the nucleonic 
current contribution while the dash-dotted curves represent the mesonic current 
contribution. The solid curves correspond to the total contribution. 



perimental efforts at COSY (see contributions by M. Wolke, D. Grzonka and 
A. Khoukaz at this meeting) are of particular interest. In any case, the value 
extracted here is compatible with the OZI value of g^NN — —(0.60 ± 0.15). 
More recent data with an improved data analysis from the DISTO collab- 
oration |3j|] exhibit a nearly isotropic angular distribution, corroborating a 
very small nucleonic current contribution. 

Fig. H shows the predicted total cross section ratio R^i^ — @pp-^pp<f>/ &pp—*p 
as a function of excess energy. For low excess energies, the predicted ratio 
is about 4 to 7 times that of the OZI estimate. At higher energies, the 
enhancement over the OZI estimate decreases to a factor of 3 or so. What 
then is the origin of this enhancement over the OZI estimate as predicted 
by our model? One source of the enhancement is in the mesonic current. 
As discussed in detail in Ref. (Tj, the violation of the OZI rule at the 4>pn 
vertex had to be introduced in order to achieve a simultaneous and con- 
sistent description of the then available data on the reactions pp — > ppui 
and pp —* pp<p. This explicit OZI violation in terms of the (fipn and uipir 
coupling constants used suggests an enhancement of around 3 in the cross 
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Fig. 4. Angular distribution of the emitted cf> meson in the cm. frame of the total 
system at a proton incident energy of 2.85 GeV. The dashed curve corresponds 
to the nucleonic current contribution while the dash-dotted curve to the mesonic 
current contribution. The solid curve corresponds to the total contribution. The 
data are from Ref . fl32l . The absolute normalization of the data has been determined 
as described in ReOTf. 



section ratio. With regard to the nucleonic current, the employed toNN and 
<f>NN coupling constants lead to results that exceed the OZI value only in 
one case: namely for the parameter set with g^N = —0.9 |7j. The corre- 
sponding enhancement factor for the cross section ratio amounts to about 2. 
It is then evident from the above consideration that the cross section ratios 
resulting from the model calculation differ significantly from those values 
implied by the employed coupling constants. Obviously, dynamical effects 
such as interferences, etc., play a rather important role here and can lead 
to a fairly large deviations from the OZI prediction within a conventional 
picture, i.e., without introducing any "exotic" mechanisms. Consequently, 
one should be very cautious in drawing direct conclusions regarding the 
strangeness content in the nucleon from such cross section ratios. The be- 
havior of the cross section ratio as the excess energy approaches zero is due 
to the finite width of the u, which prevents the oj meson production cross 
section from decreasing rapidly, as it does in the case of <f> meson. 

Evidently, the enhancement of a factor of 3 or so over the OZI estimate at 
higher excess energies in Fig. |5| is much smaller than the enhancement over 
the OZI prediction of about a factor 10 found by the DISTO collaboration 
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Fig. 5. Total cross section ratio R^/ u — o vv ^, m $l o m ^ VVUJ as a function of excess 
energy. Different curves correspond to the possible sets of parameters as determined 
in Ref.|7). The horizontal line corresponds to the OZI prediction. 

|32f at an excess energy Q w 80 MeV in <f> meson production. However, it 
is important to realize that their measurement was done at a fixed incident 
beam energy of 2.85 GeV, and therefore the corresponding excess energy of 
the produced u is already 319 MeV. Though corrections for the differences 
in the available phase space were obviously applied when extracting the 
above result, there are other effects that may influence the ratio, such as the 
energy dependence of the production amplitude, the onset of higher partial 
waves, etc., that cannot be corrected for easily. It is therefore possible 
that the actual deviation in the value of R^/^ from the OZI prediction is 
also smaller. Thus, it would be interesting to perform a measurement of the 
cross section ratio at the same (or at least similar) excess energies. From the 
theoretical side, the newer data on (ft production by the DISTO collaboration 
[j33f and, especially, the new data on u> production from COSY, which will 
become available soon p7| (see also the contribution by A. Khoukaz at 
this meeting), should already impose much more stringent constraints on 
the parameters of our model and help address better the problem of OZI 
violation and related issues in NN collisions. In fact, these new data seem 
to indicate that we overpredict the uj meson production cross section above 
Q ~ 30 MeV. This has a direct implication on the cross section ratio i2^/ w 
at higher excess energies as discussed above. 
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4.2. NN -» NNn 



We now turn to the production of pseudoscalar mesons. Here we con- 
fine our discussion to the production of n meson in NN collisions. The 
production of this meson near the threshold energy is of special interest, 
since the existing data are by far the most accurate and complete among 
those for heavy meson production and, consequently, they offer a possibility 
to investigate this reaction in much more detail than any of the other heavy 
meson production reactions. In addition to the total cross sections for the 



PP —> PPV reaction |3g, |39|, [5J [41], [4^], we have data for pn — » pnn [43 1 
and pn — > dn |42|, 44 ] . The differential cross section data for the pp — > ppr\ 
reaction jP| are available as well. Consequently, there is a large number of 
theoretical investigations on these reactions. The production of n mesons 
in NN collisions is thought to occur predominantly through the excita- 
tion (and de-excitation) of the Sii(1535) resonance, to which the 77 meson 
couples strongly. However, the excitation mechanism of this resonance is 
currently an open issue. For example, Batinic et al. ||] (see also the contri- 
bution by S. Ceci to this meeting) have found both the tt and rj exchange as 
the dominant excitation mechanism. However, they have considered only 
the pp — ► ppn reaction. Gedalin et al. j6[ and Faldt and Wilkin [l(J have 
considered both the pp — > ppn and pn — > pnn reactions. In the analysis of 
Ref.[|l0|] the pn — > dn reaction was also included. These authors || [l(| find 
the p exchange to be the dominant excitation mechanism of the 511(1535) 
resonance. In particular, it has been claimed [10] that p meson exchange 
is important for explaining the observed shape of the angular distribution 
in the pp — > ppn reaction. We also mention that, in contrast to the dom- 
inant resonance current contribution found in Refs.fi, 1(1 10 [, in a recent 



calculation of the pp — > ppn reaction by Pena et al. fly], it is found that 
the dominant contribution arises not from the 5ii(1535) resonance current, 
but from what they refer to as the short range amplitude. In our language 
this corresponds to the shorter range part of the nucleonic current. Here 
we shall report on yet another possible scenario that reproduces both the 
PP ~~ > PPV an d pn — ► pnn reactions and discuss a possibility to disentangle 
these reaction mechanisms. 

Although here we shall confine ourselves to the problem just mentioned, 
the description of n meson production in iViV collisions presents other inter- 
esting aspects. For example, the n meson interacts much more strongly with 
the nucleon than do mesons like the pion so that not only the iViV FSI, but 
also the nN FSI is likely to play an important role, thereby offering an excel- 
lent opportunity to learn about the nN interaction at low energies. In fact, 
the near-threshold energy dependence of the observed total cross section for 
n meson production differs from that of pion and 7/ production, which follow 
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the energy dependence given simply by the available phase-space together 
with the NN FSI. The enhancement of the measured cross section at small 
excess energies in 77 production compared to those in tt and rf production 
is generally attributed to the strong attractive r/N FSI. In addition to all 
these issues, the theoretical understanding of T] meson production in NN 
collisions near threshold in free space is also required for investigating the 
dynamics of the 5n(1535) resonance in the nuclear medium, possible ex- 
istence of r]NN bound states, and the possibility of using r\ to reveal the 
properties of high-density nuclear matter created in relativistic heavy-ion 
collisions. 

Our model for the r\ production current includes the nucleonic, mesonic 
and resonance currents. The mesonic current consists of the rjpp, t/wu, 
and r]a 7r exchange contributions. The resonance current consists of the 
Sii(1535), Pn(1440), and I?i3(1520) resonances excited via it, r/, p and lu 



exchange. For the NN FSI, we use the Bonn interaction [18|. For the NN 
ISI, we consider only the on-shell interaction obtained from Ref. The 



details of the calculation will be reported elsewhere [46|. The results for the 
total cross section as a function of excess energy are shown in Fig. n. The 
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Fig. 6. Total cross sections for the pp — > ppr/ (upper panel) and pn — > pnn (lower 
panel) reactions as a function of excess energy. The dashed curves correspond to the 
nucleonic current contribution while the dash-dotted curves to the mesonic current 
contribution; the dotted curves represent the resonance current contribution. The 
solid curves are the total contribution. The data are from Rcfs.[j38|, ftl| , 
[I- 
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upper panel shows the results for pp collisions while the lower panel those for 
pn collisions. The dashed curves are the nucleonic current contribution; the 
dash-dotted curves correspond to the mesonic current and the dotted curves 
to the resonance current. The solid curves are the total contribution. As 
one can see, the total cross sections are dominated by the resonance current, 
and more specifically by the strong 5n(1535) resonance (see Fig. 0). Our 
nucleonic current contributions (dashed curves) are much smaller than the 
resonance current contributions. This is in contradiction to the findings of 
Ref.||]; there, instead of the resonance current, the short range amplitude 
(shorter range part of the nucleonic current) gives a large contribution to 
the pp — > ppr] cross section. For small excess energies, our calculation 
underpredicts the data. As mentioned above, this is usually attributed 
to the r]N FSI, which is not accounted for in the present model. Note that 
the results for pn — > pnrj with excess energy Q > 50 MeV, corresponding 
to incident beam energy larger than 1.3 GeV, should be interpreted with 
caution, as no reliable NN phase shift analyses for T = states exist at 
present for energies above 1.3 GeV 



In Fig. the Sn (1535) resonance contribution to the total cross sections 
is shown (solid curves), together with the contribution from the individual 
meson exchange excitation mechanism. The dashed curves correspond to 
the 7r exchange, while the dash-dotted curves correspond to the rj exchange 
contribution. The dotted curve is due to p exchange. Although the u 
exchange is included in the calculation, its contribution is not shown here 
separately because it is much smaller than the p exchange contribution. As 
can be seen, the dominant contribution is due to the n exchange followed 
by r] exchange. The p exchange is very small. Several observations are in 
order here: 

1) The major reason for the small p exchange contribution, in contrast 
to the result of Refs.|], 10], is that in our model we have not allowed 



the vector (7^) coupling in the pNN* vertex for spin- 1/2 resonances. 
Note that for an odd-parity spin-1/2 resonance, there is an overall 
extra 7 s factor. Unlike the case of pNN vertex, the presence of this 
coupling in the pNN* vertex leads to a violation of gauge invariance, 
which is especially relevant when used in connection with the VMD. 
As mentioned in the previous section, the pNN* coupling constant in 
our model is extracted from the radiative decay, N* — » 7+ N, using the 
VMD. The simplest way of satisfying the gauge invariance constraint 
is to omit the 7^ coupling from the pNN* vertex, which we have done 
in the present work. Note that the tensor (cr^ u ) coupling is free of 
this problem. A direct consequence of omitting the vector coupling 
in the pNN* vertex is a very small p exchange contribution to the 
cross section as shown in Fie. @. Note that in Ref.g the piV5(1535) 
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Fig. 7. Same as Fig. ^, except that it shows the 5n(1535) resonance contribution 
only. The dashed curves correspond to the 7r exchange contribution while the 
dash-dotted curves to the r\ exchange contribution; the dotted curves represent the 
p exchange contribution. The solid curves show the total contribution. 



vertex is given by the 7 5 7^ coupling. An alternative to avoid the gauge 
invariance problem while keeping the 7^ term is to use a vertex of the 
form r [7^<7 2 — {tun* =p m^q^] [24|, where T~ = 7 s and T + = 1. 



In fact, Peha et al. |9[ have used a modified version of this vertex 
in conjunction with the coupling constant determined from a quark 
model [24]. They found a non-negligible contribution of the p exchange 



to the excitation of Sii(1535) in pp — > ppr]. Relevant experimental 
information should decide whether the vector coupling is required or 
not in the pNN* vertex. 

2) The 77 exchange contribution is relatively large in the present calcu- 
lation. In fact, in the case of pp — > ppr/ its contribution to the cross 
section is about half of that due to the ir exchange. The r/NN cou- 
pling strength is subject to a rather large uncertainty; the value of 
this coupling constant ranges from gNNr] ~ 2.7 to g^Nri ~ 6.4 p7| . 
The relatively large contribution of rj here result from using the r]NN 
coupling constant of g v NN = 6.14, as used in the construction of the 
Bonn NN interaction [18|. This value is close to the upper limit. 



However, the rj meson exchange in the Bonn potential [jlj| might just 
be an exchange of a (J p , T) = (0~, 0) quantum number and not of a 
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genuine rj meson. Consequently, the contribution from r\ exchange is 
subject to this uncertainty in the coupling constant. Anyway, in the 
present calculation for pp — > ppr], the rj exchange interferes construc- 
tively with the dominant ir exchange contribution, yielding the total 
contribution as shown by the solid line in Fig. |?j. For pn — ► pnr], the 
rj exchange interferes constructively with the tt exchange in the T = 1 
channel (as in the case of pp — > ppr]), but destructively in the T = 
channel due to the isospin factor —3 in the ir exchange amplitude. 

3) The correct description of both pp —> ppr] and pn — > pnr] reactions 
depends not only on different isospin factors involved for isoscalar and 
isovector mesons exchanged but, also on a delicate interplay between 
the NN FSI and ISI in each partial wave involved. While the NN FSI 
enhances the total cross section, the NN ISI has an opposite effect 
(see discussion in section II). In this connection, we mention that in 
Ref.|l(| the reduction factor due to the NN ISI is estimated to be 
about (0.77) 2 = 0.59 due to the 3 P state and (0.73) 2 = 0.53 due to 
Pi. In our calculation, however, the corresponding reduction factors 
are about 0.19 and 0.27 near the threshold energy. This large discrep- 
ancy between the results of Ref.|l(J and ours is due to the fact that, 
whereas our reduction factor is given by Eq. (|i~0|) , the reduction factor 
in Ref.|L0| is given by A = r/ 2 = ( e - /m ( 5 L)) 2 . We argue that the latter 



formula is inappropriate to estimate the effect of the NN ISI for it 
does exhibit a pathological feature: namely, when the absorption is 
maximum (t]l = 0), this formula yields A = 0, implying the total ab- 
sence of the NN elastic channel and thus not allowing the production 
reaction to occur. However, scattering theory tells us that when the 
absorption cross section is maximum, the corresponding elastic cross 
section does not vanish, but is 1/4 of the absorption cross section. 
Note that this feature is present in Eq. (|i~0|) . Furthermore, the authors 
of Ref.pl apparently have identified incorrectly the inelasticity i]l 
with cosr(p), where p is one of the two parameters (the other is the 



phase shift) given in Ref.[22|. The phase shift parametrization given 
in Ref.[2|] differs from the standard Stapp parametrization as given 
by Eq. ( |) . It is obvious that with a more appropriate estimate of the 
reduction factor A as given by Eq,(|i~0|) the result of Ref.[10| would 
underpredict considerably the cross section data. 

Fig. |8] shows the predicted angular distribution of rj in pp — > ppr] at an 
excess energy of Q = 37 MeV together with the data of Ref.^5|. Again, 
the resonance contribution (dotted curve) dominates the cross section. As 



pointed out in Ref.[10], the shape of the angular distribution of the latter 



contribution bends upwards at the forward and backward angles due to the tt 
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Fig. 8. Angular distribution of the emitted r\ meson in the cm. frame of the total 
system at an excess energy of Q — 37 MeV. The dashed curve corresponds to the 
nucleonic current contribution while the dash-dotted curve to the mesonic current 
contribution; the dotted curves represent the resonance current contribution. The 
solid curve show the total contribution. The data are from Ref . Eal . 



exchange dominance in the Sii(1535) resonance contribution. However, due 
to an interference with the nucleonic (dashed) and mesonic (dash-dotted) 
currents, the shape of the resulting angular distribution (solid curve) is 
inverted with respect to that of the resonance current contribution alone. 
As one can see, although the overall magnitude is rather well reproduced, the 
rather strong angular dependence exhibited by the data is not reproduced 
by our model. At this point one might argue that the excitation mechanism 
of the Sii(1535) resonance as given by the present model is not correct and 
that, indeed, the p meson exchange should give the dominant contribution, 
as has been claimed in Ref.[10|. However, judging the level of agreements 
between the predictions of Ref.[lC] and ours with the data, one cannot 
discard the dominance of ir and r/ exchange in favor of the p exchange 
mechanism for exciting the Sii(1535) resonance. In this connection, we 
mention that the new data from COSY which will become available soon 
shows a flat angular distribution [[47| . 

From the above considerations, we conclude that, at present, the excita- 
tion mechanism of the Sii(1535) resonance in NN collisions is still an open 
question. Indeed, we have just offered a scenario other than the p exchange 
mechanism that reproduces the available data equally well. It is therefore 
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of special interest to seek a way to disentangle these possible scenarios. In 
this connection, spin observables may potentially help resolve this issue. 



According to Ref.[lC], the p exchange contribution is expected to lead to an 



analyzing power given by 



A y = A™ ax sin(2fl) (17) 



in pp — » pprj, where A™ 11 * is positive for low excess energies, peaking at 
Q ~ 10 MeV and becoming negative for excess energies Q > 35 MeV. The 
results (dashed curves) at Q = 10 MeV (upper panel) and Q = 37 MeV 
(lower panel) are shown in Fig. |||. The corresponding predictions of the 
present model are also shown (solid curves). The different features exhib- 
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Fig. 9. Analyzing power for the reaction — > pprj as a function of emission-angle of 
r\ in the cm. frame of the total system at an excess energy of Q = 10 MeV (upper 
panel) and Q = 37 MeV (lower panel). The dashed curve corresponds to the case 
of p exchange dominance according to Ref . |l(J . The solid curve corresponds to the 
present model calculation. 

ited by the two scenarios for the excitation mechanism of the 5n(1353) is 
evident. In this connection, the COSY11 effort to measure the analyzing 
power in pp — > pprj (see the contribution by P. Winter in this meeting) is 
of great importance. We emphasize that these different results should be 
interpreted with caution. The reason for this is that, as mentioned before, 
the present model accounts for the NN ISI only in the on-shell approxima- 
tion. While this may be a reasonable approximation for calculating cross 
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sections, it may introduce rather large uncertainties in the calculated spin 
observables. 



5. Summary 

The production of heavy mesons in NN collisions has been discussed 
within the meson exchange model of hadronic interactions, paying special 
attention to the basic dynamics that determine the behavior of the cross 
sections near the threshold energy. Differences in the existing meson ex- 
change models as well as their limitations have been also discussed. Heavy 
meson production processes necessarily probe the short range dynamics, a 
domain where we have very limited knowledge so far. In this regime even 
the relevant reaction mechanisms are largely unknown. The theory of heavy 
meson production in NN collisions is still in its early stage of development. 
Successful description of these processes in terms of purely hadronic degrees 
of freedom calls for correlation of as many independent data as possible in a 
consistent way. The pp — > ppcfi reaction has been discussed as an example of 
the production of vector mesons in NN collisions. As for the pseudoscalar 
meson production, results for the n production in both the pp and pn col- 
lisions were presented. From these examples, it is clear that not only the 
pp — > ppM, but also the pn — > pnM and pn — > dM reactions should be 
investigated. Also more exclusive observables than the total cross section 
such as the spin observables should be studied. 
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